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Within the framework of General Relativity the model approach to a description of spherical 
gravitating static fluid balls with an electric charge is considered. The metric interval is written in 
Bondi's radiation coordinates. The total energy-momentum tensor as a direct sum of the perfect 
Pascal fluid energy-momentum tensor and the energy-momentum tensor of an electromagnetic 
field is chosen. The exact solution of the Einstein-Maxwell equations as an extension of similar 
solution with parabolic distribution of mass density is found. 

1. Introduction 

The problem of exact solutions finding of Einstein's equations within of the General Relativity 
does not lose a significance to this day. Special interest, despite of an exoticism, represents models 
of electric charged stars because an electrical charge is one of the few physical quantities which 
not disappears at a star blow-up and the star collapse. The basic difficulties of a calculation of 
similar models are concluded in nonlinearity of the combined equations of Einstein. It is main 
insufficiency for system description. Other difficulty is related to the solution of self-consistent of 
Einstein-Maxwell combined equations. The similar approach is realised earlier by one of authors 
in [1] where the exact solution for the charged Pascal perfect fluid with parabolic distribution 
of mass density for a neutral component of fluid is discovered, but without an introduction of a 
concrete equation of state. 

The generalisation of this model related to a modification of the mass density distribution law 
of a neutral fluid is considered in given article. Such a modification leads to a new distribution 
of the charged density into a ball model and to a new exact static spherical solution of Einstein- 
Maxwell combined equations. In this case as usually the exterior gravitational field should be 
described by Reissner-Nordstrom's solution for the charged spherical mass as the generalisation 
of exterior Swarzschild's solution. 

We make choose a static and spherical model here. The metric functions do not depend on a 
time variable. Such a model possesses spherical symmetry that removes dependence on angular 
variables. The model is considered without a rotation and a radiation. 



2. The basic mathematical expressions 

Here we write a metric in Bondi's radiation coordinates 

1 Talk presented at Proceedings of XV International Scientific Meeting "Physical Interpretations of 
Relativity Theory" (PIRT-2009), Moscow 6-9 July, 2009. 
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ds 2 = F{r)dt 2 + 2L(r)dtdr - r 2 (d6 2 + sin6 2 d(p 2 ), (1) 

where F = F(r) L — L(r) are the metric functions of a radial variable r; t is a time coordinate; 
9 and ip are angle variables; the light velocity is chosen equal to unit as a gravitational constant of 
Newton Gjq — 1. A determinant of a covariant metric tensor g a p corresponding to an expression 
(P) is det(g a/3 ) = g = —L 2 r A sin 2 8. 

We will set tetrads or an orthonormal 4-basis in a tangential space-time with using the metric 
(P) as 

9(o )M = ^; = K + \f S°; (2) 

5(2) M = -^=(Sl + isin6 Si); g (3)fl = -^7=( S l ~ 1 sin9 s lY, ( 3 ) 
9l Q) =L-'S^ g^^S^-FL-iS?; (4) 

where i is the imaginary unit and Greek indices run through 0, 1, 2, 3. 

Now we will determine basic differential 1-forms of Cartan with helping of the tetrads CE}-© 
as contractions of tetrads and the coordinate differentials 

6( Q ) = gMdaf. (6) 
We can rewrite equation ([5]) as tetrad componets 

9(0) = * pat + Ldr; = dt; (7) 

9 (2) = -j=(d9 - i sin9 dtp); 9 (3) = -j=(d0 + i sinO dip) . (8) 
Then a square of space-time interval ds 2 will be written as 

^ = s ww e (a) e(«, (9) 

where g^ a ym is a tetrad metric, which in the matrix form can be written as 



ff(a)(/3) = 9 



(0 1 

1 



Vo o -l 



(10) 



The first equations of the Cartan structure are 



J a) a e (/3) 



(11) 



where d is an exterior differential and a wedge operation A marks an exterior product. The 
equations of structure allow to find non-zero 1-forms of a connection 



(i)(o) = V^eW; W(0)(a) = (Lr)-^ 3 ); 



(12) 
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(13) 



w(i)(3) - --fWW; W(3)(2 ) - -^= C o^(e( 2 ) - e^), (w) 

a derivative with respect to radial variable r is marked by a prime. 

In general, we have for an arbitrary tetrad metric a true relationship dg/ a \/p) = U3r a \ip\ + 
%v a ), the left part of which equals zero under an investigation in our case because the tetrad 
metric components (jlOl) are constants, i.e. the 1-forms of connection has an antisymmetric 
attribute 



U(a)(p) = -W(p)(a), (15) 

This attribute allows to decrease a number of the independent 1-forms of connection. 
Further we make use of the second equations of the Cartan structure 

n (a) (fl = ^U^w a ew = d W w w + w w w A W w wi (ie) 

where fi( a )(£) = — 0( a )(,g) is a 2-form of a curvature; i2 rp)M(g) i s Riemann's tensor of the 



curvature in the tetrad notations. Non-zero components of which can be found from (fT6| and 
equals to 

F" F' L' F 2 L' 

#(1)(0)(1)(0) =~2£2 +_ 2L3" ; %)(2)(1)(2) = -4^3 ! ( 17 ) 

2' F2' F' 

«(0)(2)(0)(3) = -^3; i H0)(2)(l)(3) = - 2^2 ; ( 18 -> 

IF 

iJ (3)(2)(3)(2) = ^ - ^jp' i? (0)(3)(l)(2) = #(0)(2)(1)(3)- (19) 

According to the definition of Ricci's tensor , R( a )(p) — R ( a )(p)M = ^ R(a) {a) (p) {-/)•> 
we will write his non-zero tetrad components: 

2 2/ 

R (o)(o) = 'Ztfj;'' fl (o)(i) = 77¥ ~i?2 + ^73-772'' ( 20 ) 



FL' 


F" 


F'L' 


F' 


rL 3 ~ 


22 2 


f 2L 3 


~~ r2 2 ' 


M 1 


F 
~ Z 2 


rF' 


rF V 

V 1?T 



r 

F 2 V 

3. Einstein-Maxwell combined equations 
for charged perfect fluid 

The Einstein equations in the tetrad description with a source in the form of the energy- 
momentum tensor (EMT) are 

G (a)(p) = R (a)(P) - ^9(a)(f3)R = -xT(a)(P)i ( 22 ) 

where is the Einstein tensor; R( a )(p) is the Ricci tensor; R = R ( a j is the scalar 

curvature; x = 8ir is Einstein's gravitational constant in a chosen system of units. Resulting 
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EMT of the matter Tt a )(p) is taken as direct sum of EMT of the perfect Pascal neutral fluid 
and EMT of the electromagnetic field 



where 



J M(/3) - 1 (»)((]) + i («)(£) ' y 26 ) 



T (1)(P) = (f J ' + P) U (a)U W -P9( a )(0) = ^U {a) U m + pb( a )((3)\ (24) 



T, 



tm" = h ("^(^S + iffww f wwf ; (25) 

/i(r) is a mass-energy density; p(r) is a pressure of perfect Pascal fluid; U( Q ) = ^(q)^ — — is 
a 4-velocity in tetrad notations; fr( a )(/3) = u {a) u (&) ~ 9(a)(0) IS a 3-projector on a spacelike 
hypersurface (or 3-metric) and b^i^u^ — 0; F( a \m-\ is the tensor of electromagnetic field with 
P(a)(p) = ~F(a)(0)] a U functions here depend on a radial variable only. 
We will rewrite Einstein's equations (|22[) as 



R( a )(/3) = (T( a )(p) ~ ^9( a )(p)T^ , (26) 

using a connection between the scalar curvature and the trace of EMT, R — kT. 

Then the gravitational combined equations for a spherical symmetry can be written as four 
equations for the dimensionless radial variable x = t/Rq, which varies from zero to unit (i?o is 
the exterior ball radius) 

-^(lni)' = X T ( o) (0 ); (27) 



xL 

F 2 
2xV 



:(]nLy= X T (m)] (28) 



^(InZ)' - ^ (f» + \p - PQnLY) = - X (t (0)(1) - \ T ) ; (29) 

1 ( F xF' F . .,\ ( 1 \ 

^{- 1 + J2+J2-- XtfQnL)) = -X (T (2)(3 ) + , (30) 

d d 

where all derivatives now are taken with respect to the variable x and — = — , and a new 

or R ox 

constant is \ — — 8ttRq. 



After that we will complete our system of equations (I27l) - (|30l) with Maxwell's equations. The 
second pair of the Maxwell equations can be written in the framework of General Relativity as 

F^. v = ^=W^~gF^), v = -^f, (31) 

where j M is a density of an electric current, and a point with a comma marks a covariant 
derivation. The tensor of the electromagnetic field is set as usual over a alternation of 4-potential 
derivations F^ — A v ^ — A^^. Only one component of the tensor of the electromagnetic 
field here will be non-zero Foi = — Ao,i = — Ro A' Q . Then the Maxwell combined equations will 
be transformed to the one equation 

-l-(x 2 LF 01 y = -4TrR oJ °. (32) 

X 1j 
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Now we will take a comoving reference of frame in which the 4-velocity has the following 
components in a static case 

u m = JL= u ^ = (33) 

V^oo y/F(r) V9oo ^/F{r) 

The physical observed magnitudes will be rewritten respectively as 



Vphys = r^ w u"«" = fi(x); Pphys = p(x) = = j ^ = f^fjxj- (34) 

E»" = -F„ M «" = -Fu,^= = -^=; E phys ^ El = ^ = ^=; (35) 

V-900 y/F(x) v^oo v-r(x) 



= T#-™»u"u" = (36) 

where /Up^s is a physical observed mass-energy density; /Dpfc va is a physical observed density of 
electrical charge; Ef, hys is physical observed 3-vector of electric field strength; E p h ys = E\ is 
a radial component of physical observed of electric field strength; W e i is a physical observed of 
electric field energy density. 

Thus the Maxwell equations (|32p are rewritten in a form for the further using 



x 2 E 



4irp(x)R ^, (37) 



L J PW u Ve 

where we set a new function e = — . 

Total EMT can be written in the tetrad components as 

T (o)(o) = T (1){1) = -F(jj,+p); (38) 

T (0)(l) = 2^ ~P) + W el> T (2)(3) = r (3)(2) =P+W el . (39) 

4. Transformation of Einstein's equations 

The system of equations after a substitution of the expressions ([38|) - ((39|) into the right part 
of gravitational equations (|2"T)) - (j3"0"j) are reduced to 

£(lnL)' = f(/i+p); (40) 
x 2 

l^L)'- £ - (El + l{\nF)' -{\xxF)'{\nLA =- X (p + W el ); (41) 

-^(l-e) + £(ln^) / = -xQ(M- P ) + ^). (42) 

After that we exclude the pressure and the mass-energy density and get a linear differential 
equation with variable coefficients for function G(x) 

G" + f(x)G' + g{x)G = 0, (43) 
where G = V~F, f(x) — (In 93)', ip(x) = ^/s/x, and a coefficient g[x) equals to 
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:x (44) 



2(1 - e) + are' 2x 
2x 2 e 

An electromagnetic field influence is appeared directly through the functions s(x) and W e i(x) 
The definition of a new variable £ = £(2;), according to an expression 



(45) 

transforms (|43|) into an equation of a nonlinear spatial oscillator 

G" c ' f + f> 2 (C(:r))G = 0, (46) 
with respect to variable £. The square of "frequency" J7 2 easier to use in the form 

_ 2 d f$\ 2 Xi 



Sr = — r [-]-—W el , (47) 
ay \y J y 

because we have a difficulties with an integration of the expression (|45|) in the elementary func- 
tions. In the equation (|47|) y = x 2 , and a function $ is an analog of Newton's gravitational 
potential of an interior region of a charged fluid ball. The function $ can be found from gravi- 
tational equations through the function s as 

$ = 1 - e = | J (ji{x) + W el (x))x 2 dx = ^=J (Ml/) + W el {y))^dy. (48) 
Also an expression for pressure can be found easily from the gravitational combined equations 

XP = xWei-% + -{l-$){^F)'. (49) 

5. Boundary conditions 

The interior solution must be smoothly sewed together on the ball surface (r = Rq) with the 
exterior solution of Reissner- Nordstrom with metric ([TJ. The functions goo and <7oi are written 
in an exterior space-time as 

2m Q 2 

9oo = F R - N (r) = 1 — + 9oi = ifl-JV = 1, (50) 

where to and Q are an integral mass and an integral electric charge respectively which are 
calculated by observer on a spatial infinity. 

2m 

If we set a parameter of compactness rj = , which describes a level of gravitational 

Rq 

compression of a star then we will have the expressions for the functions F(x) = G 2 {x) and e(x) 
on a surface of ball 

F {x =i) = G\ x =x) = £(,=D = 1 - V + % = 1 - n\ (51) 

n o 

Q 2 

where r/* — rj ^ i s an effective compactness. Here we must remark that rj* < 77. 

. 

The function L automatically is continuous on a boundary r = Rq and the function $( a =i) = 

V ■ 
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The continuity requirement of a function F(x) on the ball surface leads to 

= 2(G(„i))'G ( „i) = 2rf - V- (55) 



Therefore, we have 



Now we will return to the expression for pressure (|49j) and will require of pressure's absence 
on the ball surface because there are vacuum in an exterior region. The boundary conditions 
were found for the physical functions and lead to an equation for parameters 

x W el = V - V * = (54) 
This expression is transformed in an identity, when the electric charges are absent (ry = 77*). 



6. Solution of Einstein-Maxwell equations 

Unlike paper [1], where a mass density of a neutral fluid has the parabolic distribution law 
in an interior region of the ball (Fig.l) 

fi(x) = ^o(l - bx 2 ), (55) 

we will consider the behaviour of a system for a more general mass density distribution of the 
neutral substance (Fig. 2) 

H(x) = ^{l-bx 2 f. (56) 




Figure 1. Behaviour of parabolic mass den- 
sity distribution as function of dimension- 
less radial variable for parameter b = 63/80. 



Figure 2. Behaviour of mass density distri- 
bution as function of dimensionless radial 
variable on base of parabolic mass density 
behaviour for parameter b = 63/80. 
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An expression of the electrical field energy density W e i can be generalized if we will require 
when the gravitational field disappears, i.e. Gn — > , then the expression for W e i will equal to 
the electrical field energy density for the homogeneous charged ball (p — > po = const) in a flat 
space-time 

(57 > 

where A = const. 

Therefore we will determine the electrical field energy density unlike the article [1] as 

W el = ^-\ 2 (x)x 2 , (58) 

and we will generalize also the expression of the electrical charged density so that a new expression 
of the electrical charged density with a = will be equal to the similar expression in [1] 



p = p*\Je(x) = p (l - ax 2 )\Je(x). (59) 
Then we find directly from the Maxwell equation (|37|) 



4irp Q R Q ( 3a 2 



A( S ) = (1 fx') . (60) 

We find the final expression of the electrical charged density after substitution (|60|) into (f58| 

We will find also a restriction on a parameter region a : < a < 5/3. from the expression ((54)l 
for the pressure value on the ball surface and under using equation (|6ip . However a non- negative 
behaviour of the electric charge density (fBTJf restricts else more the parameter region: < a < 1. 
We must remark that a sign of an electric charge is set with a parameter po . 

Furthermore, the boundary conditions p(x = 0) = p(0) = poy / e(fi); p(x = 1) = p(l) = 
po(l — a)y/e(l) lead us to the expression 



a = l-^hp. (62) 

If we connect parameters a and b in this case with the connection a = 80/636, then we will 
have the constancy condition of a squared "frequency" J7q = const. 

Further we can rewrite (|46l) as the equation for the harmonic spatial oscillator. 

G'l Q + Ql G = 0, (63) 

where after the using connection between parameters a and b an expression for fig can be written 
as 

n l = ^2 (170lMoa - 1760tt p 2 Q R 2 ). (64) 

A general solution can be directly written as a harmonic oscillating function for such an 
equation under condition that the squared " frequency" fig is positive 

G(C(x)) = G cos(n ■ ((x) + a ), (65) 

where «o is a phase shift. 
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The total expression of metric function goo = F = G 2 will be 



F(x) = G 2 cos 2 (fl ■ C(x) + ao), 



(66) 



and a function goi — L is found lightly from the expression L 




We can find the constants Go and ao from the boundary conditions, expressing in terms of 
known parameters 



In the summary it is necessary to remark that in the article a reduction method of the 
Einstein-Maxwell equations for a static spherical case to the equation of a nonlinear spatial 
oscillator is demonstrated. These equations were written in Bondi's radiation coordinates and 
with a source in the perfect charged Pascal fluid form. An exact solution in the framework 
General Relativity as the solution of an equation for a spatial harmonic oscillator and concrete 
relationship between the physical parameters which were included in the researching problem is 
found. Also an assumption on the mass-energy density behaviour and the electric field energy 
density in an interior region of the charged fluid gravitational ball was introduced. This solution 
is a generalization of an earlier finding solution in [1]. 
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(67) 



tan(fi C( x = 1) + ao) = - 



2?7* - r\ 



(68) 
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